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ABSTRACT 

Perhaps as many as 30 parallax microlensing events are known, thanks to the efforts of the 
MACHO, OGLE, EROS and MOA experiments monitoring the bulge. Using Galactic models, 
we construct mock catalogues of microlensing light curves towards the bulge, allowing for the 
uneven sampling and observational error bars of the OGLE-II experiment. As a working def- 
inition of a parallax event, we require the improvement Ay 1 on incorporating parallax effects 
in the microlensing light curve to exceed 50. This enables us to carry out a fair comparison 
between our theoretical predictions and observations. The fraction of parallax events in the 
OGLE-II database is around ~ 1 per cent, though higher fractions are reported by some other 
surveys. This is in accord with expectations from standard Galactic models. The fraction of 
parallax events depends strongly on the Einstein crossing time fa, being less than 5 per cent 
at ?e ~ 50 days but rising to 50 per cent at > 1 yr. We find that the existence of paral- 
lax signatures is essentially controlled by the acceleration of the observer normalised to the 
projected Einstein radius on the observer plane divided by f|. The properties of the parallax 
events - time-scales, projected velocities, source and lens locations - in our mock catalogues 
are analysed. Typically, ~ 38 per cent of parallax events are caused by a disk star microlensing 
a bulge source, while ~ 33 per cent are caused by a disk star microlensing a disk source (of 
these disk sources, one sixth are at a distance of 5 kpc or less). There is a significant shift 
in mean time-scale from 32 d for all events to ~ 130 d for our parallax events. There are 
corresponding shifts for other parameters, such as the lens-source velocity projected onto the 
observer plane (~ 1110 km s _I for all events versus ~ 80 km s for parallax events) and 
the lens distance (6.7 kpc versus 3.7 kpc). We also assess the performance of parallax mass 
estimators and investigate whether our mock catalogue can reproduce events with features 
similar to a number of conjectured 'black hole' lens candidates. 

Key words: gravitational lensing - Galaxy: bulge - Galaxy: centre - Galaxy: kinematics and 
dynamics. 



1 INTRODUCTION 

Thousands of microlensing events in the Local Group have been 
discovered by various collaborations, such as MACHO (e.g., Al- 
cock et al. 2000), OGLE (e.g., Wozniak et al. 2001, Udalski 2003), 
MOA (e.g., Sumi et al. 2003), EROS (e.g., Afonso et al. 2003) and 
POINT-AGAPE (e.g., Paulin-Henriksson et al. 2003). The vast ma- 
jority of these are toward the Galactic centre, many of which were 
discovered in real-time. 1 The database of microlensing events pro- 
vides a unique mass-selected sample to probe the mass function 
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1 For example, see http://www.astrouw.edu.pl/~ogle/ogle3/ews/ews.html 
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of lenses and the mass distribution and dynamics of the Galaxy 
(see Paczynski 1996 and Evans 2003 for reviews). Unfortunately, 
the mass of the lens can not be unambiguously determined for 
most microlensing events because of degeneracies. However, for 
the so-called exotic microlensing events, which include the finite 
source size events (Witt & Mao 1994; Gould 1994) and parallax 
microlensing events (Gould 1992), the degeneracies are partly or 
wholly broken. Astrometric microlensing offers another exciting 
possibility to determine the lens mass (e.g., Walker 1995). As- 
trometric measurements may become feasible with the VLT inter- 
ferometer (Delplancke, Gorski & Richichi 2001) and future satel- 
lite missions such as the Space lnterferometry Mission (Paczynski 
1998) and GAIA (Belokurov & Evans 2002). 

The standard microlensing light curve follows a characteristic 
symmetric curve (e.g., Paczynski 1986). However, this is based on 
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the assumption that the relative motions among the observer, lens 
and source are all uniform and linear. This assumption is clearly 
wrong in principle as we know the Earth revolves around the Sun, 
and furthermore, the lens and source may be in binary systems 
of their own. However, as most microlensing events last of the 
order of weeks, the effects of acceleration are not noticeable in 
most events. Nevertheless, for some events, the resulting departures 
from the standard curve are clearly visible. They can range from a 
slight asymmetry to dramatic multiple peak behaviour (Smith et 
al. 2002a). The significance of these so-called parallax events is 
that they allow an additional constraint to be placed on the lens 
mass. The parallax events are biased towards (i) more massive, (ii) 
slow-moving or (iii) closer lenses. The first of these biases means 
that parallax events offer a powerful way to detect stellar remnants, 
such as neutron stars and stellar mass black holes (Agol et al. 2003). 
Combined with other exotic effects, such as finite source size effect, 
one can derive the lens mass uniquely (e.g., Jiang et al. 2005). 

A systematic survey of parallax events of the MACHO 
database was performed by Becker (2000). Bennett et al. (2002a) 
subsequently published the most convincing long-duration events 
from this survey. Smith, Mao & Wozniak (2002b) have searched 
systematically for the parallax events in the 3-yr OGLE-II database. 
Parallax events have also been found serendipitously in the MOA 
(Bond et al. 2001) and EROS databases (Afonso et al. 2003). Table 
Qis a compendium of the good and marginal parallax candidates in 
the direction of the Galactic bulge. The fraction of parallax events 
in the microlensing database ranges from around ~ 1—10 per cent, 
although some events are more convincing than others. Additional 
parallax events have been identified by Popowski et al. (2004), al- 
though no model fits were presented for these candidates. Important 
advances not listed in this table include the detection of the parallax 
effect in the binary lens events EROS BLG-2000-5 (An et al. 2002) 
and OGLE-2002-BLG-069 (Kubas et al. 2005), and in the Large 
Magellanic Cloud single-lens event MACHO LMC-5 (Gould, Ben- 
nett & Alves 2004). For all of these events, the parallax effect leads 
to an accurate determination of the mass of the lens, when com- 
bined with the other astrometric or photometric data. 

The question naturally arises whether the observed fraction is 
consistent with theoretical expectations. This is an important ques- 
tion as the number of parallax events depend on the mass function 
and kinematics of lenses. There have been several previous studies 
(Buchalter & Kamionkowski 1997; Bennett et al. 2002a). However, 
these studies have some deficiencies. For example, they all used 
regular samplings and uniform simulated errors. Even sampling is 
clearly a gross simplification as there are significant gaps in the 
observational data, in particular, the annual period (~ late-October 
- mid-February) during which the bulge cannot be observed. This 
gap can be especially important for events of sufficiently long du- 
ration, since the asymmetric nature of the parallax signal will be 
more difficult to detect if a significant part of the event lies within 
this gap. 

The purpose of this paper is to conduct a study of parallax 
events using Monte Carlo simulations for the OGLE-II experiment. 
We will explicitly account for the uneven sampling, simulate realis- 
tic error bars and adopt the same event selection criteria as Wozniak 
et al. (2001), in order to make a much better comparison between 
observations and theoretical predictions. The outline of the paper 
is as follows. We first present the details of our simulations in Sec- 
tion [2] We then make mock catalogues of microlensing events to- 
wards the Galactic bulge and compute the fractions of observable 
parallax events in Section We analyse the properties of the par- 
allax events in our mock catalogues in Section [4] This section in- 



cludes a number of subsections dealing with event parameters <4. II 
- 14.31 . mass estimators I4.4i and long duration black-hole can- 
didate events I4.5i . We finish with a summary and discussion of 
strategies for future parallax surveys in Section|5| 



2 MONTE CARLO SIMULATIONS 
2.1 Parallax Microlensing Events 

The observable quantities for parallax events are the Einstein radius 
crossing time t E and the projected velocity v on the observer plane. 
The former quantity is, 



t - rE 

Vi 



4GM 



-D s .x(l - x), 



(1) 



where r E is the Einstein radius, M is the mass of the lens, v ± is the 
speed of the lens transverse to the observer-source line of sight, £>s 
is the source distance, Z\ is the lens distance, while x = Z\/Ds is 
the ratio of the distance of the lens to the source. The latter quantity 
is (Gould 1992), 
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where v L , v s and vq are the velocities of the lens, source and the 
Sun transverse to the line of sight. The Einstein radius crossing 
time f E is always measurable for well-sampled microlensing events; 
however, the projected velocity is only measurable if the magnifi- 
cation fluctuations caused by the motion of the Earth are substan- 
tial. The observables (f E , v) can be used to construct another useful 
quantity, namely the Einstein radius projected on the observer plane 



h = vt E . 



2.2 The Galactic Model 



(3) 



We assume that the sources and lenses may lie either in the Galac- 
tic disk or the bulge. In practice, there may be some contamina- 
tion from sources in the Sagittarius dwarf galaxy (e.g., Evans 1995; 
Cseresnjes & Alard 2001). Bennett et al. (2002a) have argued that 
this may be particularly important for the long duration events. 
Nonetheless, the structure of the disrupting Sagittarius dwarf is too 
irregular and uncertain for reliable modelling and so we do not in- 
clude it as a source population. In any case, Cseresnjes & Alard 
argue that the contribution from Sagittarius source events is less 
than 1 per cent for the OGLE-II fields. 

We adopt a value of 8.5 kpc for the Galactocentric distance of 
the Sun. Using standard cylindrical polars (R, z), the density law of 
the Galactic disk is (e.g., Binney & Evans 2001) 



Pl(R, z) = p exp - 



R 



0.435 sech* 



3 kpc 

z 



+ 0.565 exp - 



kl 



(4) 



(220 pc/ ~ r \ 440 pc J 

This uses the vertical profile found by Zheng et al. (2001) in 
their studies of disk M dwarfs with the Hubble Space Telescope. 
By counting stars within 5 pc of the Sun (which can be detected 
through their large proper motions) and using Hipparcos paral- 
laxes JahreiB & Wielen (1997) find that stars contribute 3.9 x 10~ 2 
Mq pc~ 3 to the mass density at the plane, which sets the local mass 
density and hence the overall normalisation po. This is an accurate 
representation of the local disk, embodying information from local 
star counts and stellar kinematics. The velocity distribution of disk 
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Table 1. Data for the known parallax events towards the Galactic bulge. Two events (OGLE-1999-BUL-32/MACHO-99-BLG-22 and MACHO-96-BLG- 
12/EROS-BLG-12) have two separate sets of best-fit parameters since they have been observed and modelled by two different collaborations. Errors are not 
given, but can be found in the relevant reference. Values in parenthesis indicate the presence of two degenerate parallax fits. Additional MACHO events 
have been presented in Becker (2000), although the more convincing events from this sample are presented below, since they also appeared in Bennett et al. 
(2002a,b). 



Lightcurve (d) v (km s ) (au) Reference 



OGLE sc6.2563 


71.58 


99.3 


4.1 


Smith et al. (2002a) 


OGLE sc20.5748 


78.2 


53.9 


2.4 


Smith et al. (2002a) 


OGLE sc27.3078 


124 


23.5 


1.7 


Smith et al. (2002a) 


OGLE sc33.4505 


194 


56.9 


6.4 


Smith et al. (2002a) 


OGLE sc4 1.3299 


98 


41.1 


2.3 


Smith et al. (2002a) 


(J(jLb sc43_836 


45.1 


58.0 


1.5 


Smith et al. (2002a) 


OGLE sc26.2218 a - b 


39.53 


167 


3.8 


Smith, Mao, Wozniak (2003a) 


OGLE-1999-BUL-19 


372.0 


12.5 


2.7 


Smith et al. (2002b) 


OGLE-1999-BUL-32 


640 


79 


29.1 


Mao et al. (2002) 


- a.k.a. MACHO-99-BLG-22 


560 


75 


24.3 


Bennett et al. (2002b) 


OGLE-2000-BUL-43 C 


156.4(158.2) 


40.1 (52.4) 


3.6 (4.8) 


Soszyhski et al. (2002) 


OGLE-2003-BLG-238 a 


38.2 


652.7 


14.4 


Jiang et al. (2004) 


OGLE-2003-BLG-175 d 


-63 


-141 (-106) 


-5.1 (-3.9) 


Ghosh et al. (2004) 


MACHO- 104-C 


110 


77 


4.9 


Bennett et al. (2002a) 


MACHO-96-BLG-5 


485 


30.9 


8.7 


Bennett et al. (2002a) 


MACHO-96-BLG-12 


147 


47.5 


4.0 


Bennett et al. (2002a) 


-a.k.a. EROS-BLG-12 


145.6 


43.7 


3.7 


Afonso et al. (2003) 


MACHO-98-BLG-6 


245 


79 


11.2 


Bennett et al. (2002a) 


MACHO-99-BLG-1 


115.5 


43.9 


2.9 


Bennett et al. (2002a) 


MACHO-99-BLG-8 


120 


62 


4.3 


Bennett et al. (2002a) 


EROS-BLG-29 


108.3 


69.8 


4.4 


Afonso et al. (2003) 


MOA-2000-BLG-11 


69.7 


42.5 


1.7 


Bondetal. (2001) 


MOA-2003-BLG-37 d 


-43 


-70 (-50) 


-1.7 (-1.3) 


Parketal. (2004) 



"This event also exhibits finite source signatures. 

b The parallax detection becomes marginal if the possibility of small, negative blending is allowed in the fitting. 

c Analysis of additional EROS data for this event has been able to discriminate between the two degenerate sets of parameters, showing that the fit with 
t E = 158.2 d is unfeasible (Le Guillou 2003). 

d Only approximate values are quoted as the parallax parameters have been found to suffer from additional degeneracies. 



stars is taken as a Gaussian with a mean <v> = (0,214,0) km s -1 
in cylindrical polar coordinates. The disk velocity dispersion is di- 
agonalised along the same cylindrical polar coordinate axes with 
cr RR = 21 km s~\ cr^ = 34 km s -1 and av. = 18 km s _1 (Edvards- 
son et al. 1993). 

The density law for the bulge deflectors is 



PL(x',y',z) = 
Po exp ] - 



990 pc 



385 pc 



250 pc 



1/2 



(5) 



with the major axis (;t'-axis) in the Galactic plane and oriented at 
~24° to the line of sight. This is the E2 model, as suggested by 
Dwek et al. (1995) and subsequently modified by Stanek (1997). It 
is a good fit to the near-infrared photometry of the bulge as seen 
by the COBE. The normalisation p is set to be 4.46 Mq pc~ 3 , 
which is obtained by setting the total mass within 2.5 kpc of the 
Galactic centre to be 1.5 x 10 10 Mq. The velocity distribution of the 
bulge stars is a Gaussian about zero mean. The velocity dispersion 
tensor is diagonal in the Cartesian coordinates along the axes of 
the triaxial bulge with <Tyy = 114 km s -1 , 0y T < = 86 km s -1 and 
<r a = 70 km s" 1 (c.f., Han & Gould 1995; Evans & Belokurov 
2000). 

The mass function if>(M) - or the number of stars per unit mass 
- is the multi-part power-law taken from equation (4) of Kroupa 
(2002). This is an initial mass function and so coincides with the 
present-day mass function below w 1M Q . Above 1M Q , the index of 



the power-law is adjusted to give a good representation of the lumi- 
nosity function in the OGLE-II fields; this gives a steep power-law 
of -7. The luminosity function corresponding to the mass func- 
tion is computed in the following way. Baraffe et al. (1998) pro- 
vide mass-to-light coefficients in different bands for stars of differ- 
ent ages with different chemical abundances. We assume a stellar 
age of 5 Gyr and solar metallicity, then compute the transformation 
from mass function to luminosity function by numerical differentia- 
tion of this data. In our simulations, masses are generated between 
0.05M o and 5M Q . We do not generate masses above 5M Q , as the 
probability distribution is a sharply decreasing function of mass. 
All the stars with mass below 0.08M Q are treated as dark. 



2.3 Simulation Algorithm 

To simulate microlensing events, we must pick the event parame- 
ters from the probability distribution 



P(t, b, £> s , m) oc D 2 S p s (£, b, D s )f(m, C, b, D S )T(£, b, D s ) 



(6) 



Here, p s is the density of sources, *F is the luminosity function and 
T is the microlensing rate at the source location (£,b,Ds). We as- 
cribe the source population to the bulge or disk according to the 
density at this (£, b, D s ) and choose random velocity components 
according to the source population. We note that the source veloc- 
ity probability distributions are separable. 

The luminosity function *F depends on position because of ex- 
tinction. The V-band extinction at any location is calculated using 
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Figure 1. Four sample light curves from our mock catalogues. These light curves illustrate the different cuts that are applied to the catalogue to differentiate 
between no-, weak-, moderate-, and strong-parallax signatures (top left, top right, bottom left and bottom right, respectively). By our definitions, a no-parallax 
event is one with A^ 2 < 10. The solid line denotes the true light curve, while the dotted and dashed lines correspond to the best-fitting standard and parallax 
models, respectively. For each light curve we also plot the residuals from the true light curve (top panels). The error-bars have been omitted to aid clarity, 
although a typical error bar has been included - note that this is not an actual data point from the event. 



Drimmel & Spergel's (2001) extinction law and translated into /- 
band extinction using A V /E(V - I) = 2.1 (e.g., Popowski 2001). 
We generate sources within the magnitude range 13.6 < I < 21.0. 
The lens mass is generated from 

P(M) oc M 1/2 0(M). (7) 

The flux contributed by the lens is calculated and the microlens- 
ing event is retained only if / < 19. The cut for objects fainter 
than / = 19 is applied to reduce the problem of blending by faint 



background sources. Our algorithm therefore only takes the blend- 
ing by the lens light into account and applies a cut in magnitude to 
minimise the effects of blending by faint background sources. Al- 
though, in practise, the flux from the source can be contaminated by 
light from other nearby stars, in Section|5|we show that our results 
are effectively unchanged if we incorporate a simple distribution 
for this additional blended light. 

The differential rate can be written out explicitly as 
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p L (£&,D L )F(v i )2v J .r E (M,A.) 



(8) 



d 4 r _ <f>(M) 
d 2 v L dD L dM " ~~(M) 

where p L is the total density of lenses, F(v ± ) is the distribution of 
relative transverse lens velocities, (M) is the average mass while 
v ± is the relative transverse velocity of the lens, as defined in equa- 
tion J2j. So, for the lenses, we choose the distance Dl and velocity 
v L from the probability distribution: 



P(D L ,v L )oc ^D l (D s -D l )Pl(D l )F(v ± )v 



(9) 



The lens population is chosen in this step by comparing the relative 
densities of disk and bulge lenses at this location. Finally, an impact 
parameter (defined as the minimum separation between the lens 
and the Sun-source line of sight, in units of the Einstein radius) 
is generated in the range 6 [-2,2]. This approach to selecting the 
impact parameter means than a fraction of events will be incorrectly 
omitted for r E < 1, since the peak magnification depends on the 
separation between the lens and the Earth-source line of sight, not 
the Sun-source line of sight; however, this problem is negligible 
since only 0.3 per cent of events from our model have f E < 1. 

This prescription gives us a microlensing event for our mock 
catalogue. 




2.4 Mock Light Curves and Catalogues 

The characteristics of the light curves were chosen to match the 
520 event OGLE-II Difference Image Analysis (DIA) microlens- 
ing catalogue of Wozniak et al. (2001; see also Wozniak (2000) for 
a detailed description of the OGLE-II data). These OGLE-II cata- 
logues are constructed from three years of Galactic bulge observa- 
tions, with observations taken in the /-band once every few nights 
during the bulge season (typically mid-February until the end of 
October) resulting in between 200 and 300 observations per light 
curve. The OGLE-II experiment consists of 49 bulge fields cover- 
ing a total area of approximately 10 deg 2 (see fig. 1 of Wozniak 
et al. 2001). Two of the fields are monitored much less frequently 
and a further three fields had no observations during the first sea- 
son. The limiting magnitude of the experiment is / a 20 and the 
saturation limit is / » 11.5. 

To obtain the time sequence of observations we randomly se- 
lect 100 OGLE-II light curves for each field from the variable star 
light curves of Wozniak et al. (2002), and then pick one of these 
time-series at random for each simulated event in the field. We 
then calculate the flux for each epoch. To do this, we assume that 
the photometric errors are Gaussian and scale according to the fol- 
lowing empirical relation derived from the variable light curves of 
Wozniak et al. (2002), 

o- F = 5.04 + 3.58 x 10~ 3 F lM . (10) 

Here, the flux (F) and the error on the flux (<r F ) are in units of 
10 ADU. These fluxes can be converted into /-band magnitudes 
through, 

/(?) = 23.35 - log F(t). (11) 

As stated in Section l2~3l we simplify our simulations by incorpo- 
rating the blended flux from the lensed star only, i.e. it is assumed 
that all of the observed flux comes from either the lensed source or 
the lens and not from any other nearby stars. 

As a final step, we implement the microlensing detection cri- 
teria of Wozniak et al. (2001). These criteria were employed to 
discriminate against variable stars in the OGLE-II DIA catalogue 
and must be applied to our mock catalogues in order to select only 



Figure 2. Comparison between our mock catalogue (solid line) and the ob- 
served OGLE-II catalogue of 33 bright events (data points) from Sumi et 
al. (2005), where we have assumed Poisson errors. The top panel shows 
the time-scale distribution and the lower panel shows the amplitude of the 
magnification at the peak of the event. To make a fair comparison with the 
catalogue of Sumi et al (2005) we have only plotted events from our mock 
catalogue that are unblended, have magnifications greater than 3/V5 and 
have baseline magnitude brighter than / = 17. 



S/Ncut 


Percentage of parallax events 




Weak parallax 


Moderate parallax Strong parallax 




(A X 2 > 10) 


(AX 2 > 50) (Ax 2 > 100) 


5 


4.6 


1.4 0.8 


15 


7.1 


2.3 1.5 


30 


10.0 


3.6 2.3 



Table 2. The percentage of parallax events using the criteria defined in Sec- 
tion l3.1l as a function of S/N cut. 



events with noticeable brightening and a sufficiently constant base- 
line. 

To assess the credibility of our mock catalogue, we compare it 
to the observed OGLE-II catalogue of 33 bright events from Sumi 
et al. (2005). This catalogue of Sumi et al. (2005) is based on events 
that lie within an extended Red Clump Giant region and are se- 
lected to be unblended. We choose this catalogue rather than the 
520 event catalogue of Wozniak et al. (2001) as no detailed mod- 
elling was carried out on this larger sample. In addition, the cata- 
logue of Wozniak et al. (2001) is affected by blending from unre- 
lated stars near to the source, whereas our mock catalogues only 
consider blending from the lens. In Fig.|2|we compare two proper- 
ties, the r E distribution and the distribution of amplification at the 
peak. As can be seen from this figure, the two catalogues are in 
good agreement. For example, the values of the mean t E are 29.4 
d and 28.9 d for the mock catalogue and the observed catalogue, 
respectively. 
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Figure 3. The distribution of Ax for our simulated catalogue (solid) and 
for the same catalogue but with no parallax signatures (dashed). All events 
have S/N > 5. 

3 FREQUENCY OF PARALLAX EVENTS 

3.1 Working Definition of A Paralfax Event 

As a first step in analysing our mock catalogues, we fit all events 
with both the five-parameter blended Paczynski light curve and 
seven-parameter blended parallax model. The improvement af- 
forded by the parallax model is recorded as the improvement in 
X 2 ■ We also calculate a measure of the signal-to-noise ratio (S/N) 
for each light curve, defined as the minimum value of F/cr F for the 
three points bracketing the maximum amplification. 

Before comparing theoretical predictions with observations, 
we must first investigate what level of improvement in x 2 is re- 
quired to classify an event as a parallax event. This is important 
because it is conceivable that problems such as scatter in the data 
could be misidentified as parallax signatures. 

To test this, we construct a catalogue of events using our 
Galactic model but generate light curves with no parallax signa- 
tures. We then fit these events with both the standard and paral- 
lax models and investigate the distribution of Ax 2 = xtim ~ Xpam- 
Fig. shows the distribution of Ax 2 for both the parallax and no- 
parallax catalogues. The no-parallax distribution shows a rapid de- 
cline in the fraction of events with Ax 2 > 10 (less than 1 per cent 
have Ax 2 > 10). This rapid decline indicates the minimum cut that 
should be applied in order to identify parallax events. However, 
when dealing with real data, we must also bear in mind that there 
may be contamination from other effects, such as binary signatures 
and/or problems with the data. 

Throughout this paper, we adopt three different designations 
for parallax events: convincing or strong events, with Ax 2 > 100; 
moderate events, with Ax 2 > 50; and marginal or weak events, 
with Ax 2 > 10. (If one uses the F-test for the significance of 
parameters [see Smith et al. 2002b], these limits correspond to 
log p F < -2,-10,-20, respectively). Unless otherwise stated, 
when we refer to parallax events, the 'moderate' criterion is im- 
plied, i.e. Ax 2 > 50. fn addition, unless otherwise stated, we restrict 



Strong Parallax 



Moderate Parallax 




0.2 0.4 0.6 0.8 



^ r E I ( r E,flt r E.truc)/ r E.truo 

Figure 4. The distribution of fractional error in the recovered projected Ein- 
stein radius, i.e. <5re = |{?Wflt - ' s E,tme)/^E,trucl- This is shown for our three 
levels of parallax events, strong (long dash), moderate (solid line) and weak 
(short dash). The large fraction of weak parallax events with S?e ~ 1 are 
caused by events with fE,m ~ 0. 

ourselves to good quality parallax events with S/N > 5 (as defined 
at the beginning of this section). In Fig.Q we show four sample 
light curves from our mock catalogues. These light curves show 
clearly how the standard model becomes increasingly less able to 
fit the mock data as Ax 2 becomes greater. They also illustrate that 
gaps in the data (in particular, the three month gap between ob- 
serving seasons) significantly affect our ability to identify parallax 
events. 

The issue of classifying parallax events is further complicated 
by the fact that even though an event may be displaying parallax 
signatures, this does not necessarily imply that the parallax param- 
eters can be recovered with a high degree of accuracy due to the 
degeneracies that are inherent in the parallax formalism (see Gould 
2004 and references therein for details of the various types of de- 
generacies - both continuous and discrete - that can affect parallax 
events). We illustrate this in Fig. |4]by investigating the distribu- 
tion of fractional error in recovered projected Einstein radius, i.e. 
8tr = K/E.m - ' ; E,true)/' : E,triiel- We find that the fraction of events with 
8r£ < 0.3 is 64, 53 and 23 per cent for our strong, moderate and 
weak parallax events, respectively. The large fraction of weak par- 
allax events with STe ~ 1 are caused by events with r^g, ~ 0. A 
comprehensive investigation into the nature of the errors for the fit- 
ted parameter f E is beyond the scope of this paper. 

3.2 Fraction of Parallax Events 

We now analyse the mock catalogue in order to calculate predicted 
fractions of parallax events. To do this, we compute the total num- 
ber of events that have Ax 2 greater than a given value. These frac- 
tions are plotted in Fig.|5|as a function of the Ax 2 cut for two differ- 
ent choices of S/N cut. This figure shows that for moderate parallax 
events (i.e. Ax 2 > 50) our catalogues predict a total of ~ 1 .4 per cent 
from all events with S/N > 5 or -3.6 per cent for S/N > 30. We 
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Figure 5. The cumulative fraction of parallax events in the mock catalogue, 
i.e. the fraction of parallax events for different cuts of A% 2 . The thick and 
thin lines are for events with S/N > 5 and S/N > 30 respectively. 



Figure 6. The lower panel shows the distribution of event time-scales for all 
events (dotted line) and parallax events (solid line) in our mock catalogue. 
The upper panel shows the corresponding fraction of parallax events as a 
function of time-scale for strong, moderate and weak parallax events (long 
dash, solid and short dash, respectively). Note that these plots use the true 
event time-scale, not the fitted event time-scale. 



tabulate our findings in Table [2] For Axl ul > 10, the slope of these 
relationships can be approximated by the following power-laws, 

P(A X 2 > AxL) = 0.22(4^,,)-° 71 for (A^ ul > 10, S/N > 5) (12) 
P(A X 2 > Axlj = 0.44(A X 2 U1 )- 064 for (A X 2 cut > 10, S/N > 30). (13) 

We have also computed the fraction of parallax events on 
replacing our standard Galactic disk model with a maximal disk 
model, in which the local disk surface density is increased to 100 
M pc~ 2 . The rationale for this is that local disk stars are good can- 
didates as lenses for parallax events, and so the fraction may be 
a good diagnostic of the local disk density. However, while this 
change does produce slightly more parallax events, the effect is 
small and dwarfed by the observational uncertainties. 

The total fractions from our mock catalogues are in broad 
agreement with a previous study by Buchalter & Kamionkowski 
(1997), which estimated that ~ 1 per cent of microlensing events to- 
wards the bulge should exhibit noticeable parallax signatures. The 
results from Fig. [5] are also in rough agreement with the fraction 
of parallax events from different microlensing collaborations, such 
as MACHO, OGLE and MOA. From a catalogue of 321 events 
from a 7-yr survey by the MACHO collaboration, approximately 
~2 per cent exhibited convincing parallax signatures and approx- 
imately ~3 per cent exhibited weak parallax signatures (Bennett 
et al. 2002a). The MOA collaboration found one parallax events 
(Bond et al. 2001) in their 20-event catalogue from bulge obser- 
vations in the year 2000, i.e. ~5 per cent. The fraction of paral- 
lax events in the OGLE-II catalogues can be found in two studies: 
Smith et al. (2002b) found 1 convincing parallax event from a sam- 
ple of 512 events, although a number of parallax events were sub- 
sequently found to have been omitted from this 512 event sample 
(e.g. Mao et al. 2002; Smith et al. 2002a); and Sumi et al. (2005) 



found 1 parallax event 2 in their sample of ~ 30 bright red clump gi- 
ant sources. Although the above results seem to be consistent, a re- 
cent study by the EROS collaboration (Afonso et al. 2003) reported 
a much higher fraction of parallax events; from a total of 16 red 
clump giant events they identified 2 convincing parallax events, i.e. 
12.5 per cent. We return to the issue of these EROS parallax events 
in the discussion. However, the level of discrepancy (or agreement) 
between these fractions must be taken with caution as all are sub- 
ject to small number statistics. In addition, such comparisons are 
affected by the different properties of the catalogues, i.e. the ob- 
serving strategy and duration of each project, etc. When comparing 
such fractions, we should also bear in mind the fact that our mock 
catalogues are restricted to events with /-band baseline magnitude 
brighter than 19 mag. 



4 THE PROPERTIES OF THE PARALLAX EVENTS 

In this section, we investigate various properties of parallax events 
in our mock catalogue in order to understand their nature. We utilise 
the criteria defined in Section l3~T1 to refer to strong, moderate and 
weak parallax events and enforce a minimum S/N of 5. We deal 
with the event timescales, lens and source distances, and projected 
velocities in turn. Table [3] summarises the mean and median val- 
ues for these quantities. We also assess the performance of parallax 
mass estimators and investigate whether our mock catalogue can 
reproduce events with features similar to a number of conjectured 
'black hole' candidates. 



2 This parallax event is the same as the one discovered in Smith et al. 
2002b. 
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Mean (median) values 





All events 


Weak parallax 


Moderate parallax 


Strong parallax 


te(d) 


32.1 (22.9) 


91.2 (73.4) 


130.4(105.8) 


141.6(112.0) 


h. (au) 


10.3 (8.3) 


6.4 (5.2) 


4.8 (3.9) 


4.5 (3.8) 


v (km s -1 ) 


1108.6 (609.0) 


247.2(113.9) 


80.1 (61.4) 


69.6 (52.8) 


£>s (kpc) 


9.7 (9.2) 


9.6 (9.0) 


9.2 (8.6) 


9.1 (8.6) 


D L (kpc) 


6.7 (7.0) 


4.7 (4.4) 


3.7 (3.1) 


3.5 (2.9) 



Table 3. Mean (median) values for the microlensing parameters discussed in SectionRI 




Figure 7. The upper panels show the fraction of events that are classified as parallax as a function of lens/source distance for strong, moderate and weak 
parallax events (long dash, solid and short dash, respectively). The lower panels show the distributions of lens and source distances for all events (dotted line) 
and parallax events (solid line) in our mock catalogue. The vertical dashed lines at 8.5 kpc denote the centre of the Galaxy in our simulations. 



4.1 Event Time-scale 

Fig.[6]shows the distribution of the true event time-scale 3 (i.e. not 
the fitted time-scale). As expected, the parallax events all lie in the 
long duration tail of the time-scale distribution. While the t E distri- 
bution for all events in our mock catalogue peaks around t E ~ 20 d, 
parallax events peak around t E « 100 d. Although we might expect 
parallax events to be detected at very large time-scales, few are 
seen with t E > 500 d. This is because our simulated light curves are 
based on the 3-yr dataset from OGLE-II. The microlensing event 
detection drops off sharply for t E > 500 d, as such long events do 
not pass the constant baseline criterion necessary to differentiate 
microlensing events from variable stars (Wozniak et al. 2001). 



Given the quality and sampling of the OGLE-II data, it is 
unlikely that we can detect parallax signatures for events with 
t E £ 50 d, even for weak parallax events. The upper panel of Fig.|6| 
shows that at t E ss 50 d, the fraction of events that display parallax 
signatures is less than ~ 5 per cent (0.7, 1.5, and 7.2 per cent for 
strong, moderate and weak parallax events). The prospect of ob- 
taining or constraining the parallax parameters for such short du- 
ration events can be helped by the more-frequent sampling and/or 
improved photometric accuracy, such as that afforded by the nu- 
merous microlensing follow-up networks (e.g. Jiang et al. 2005); 
however, it is known that such short duration events are affected by 
significant degeneracies (e.g. Gould 1998). 



3 It should be noted that the underlying distribution of time-scale is mainly 
controlled by the choice of mass function (see, for example, the left hand 
panel of figure 10 in Alcock et al. 2000); however, the purpose of our work 
is to investigate parallax events and so we do not concern ourselves with 
undertaking a detailed analysis of the overall form of the time-scale distri- 
bution. 



The fraction of events that display moderate parallax signa- 
tures exceeds 50 per cent for events with t E <; 1 yr. However, the 
fact that the fraction of parallax events only reaches ~50 per cent 
for t E « 1 yr shows that even for such long duration events, the 
presence of strong parallax signatures is not guaranteed. 
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Figure 8. The lower panels show the distribution of the Einstein radius projected into the observer plane left) and the projected lens velocity on the observer 
plane (v; right) for all events (dotted line) and parallax events (solid line). The upper panels show the fraction of events that display parallax signatures as a 
function of r E (left) and v (right) for strong, moderate and weak parallax events (long dash, solid, and short dash, respectively). Note that these plots use the 
true values for r E and v, not the fitted values. 



4.2 Lens and Source Distances 

Fig. shows the distributions of lens and source distances. The 
lensed sources are preferentially on the far side of the Galactic cen- 
tre (see also Table [3] which shows that the mean source distance 
for all events is 9.7 kpc), as these sources have more foreground 
lensing stars, i.e. the optical depth is greater for sources behind the 
Galactic centre. They will appear fainter (Stanek 1995) and may 
have different proper motions from the overall population of stars 
(Mao & Paczynski 2002). 

From Fig. [7] it can be seen that nearby lenses are favoured 
for parallax events. This is due to the fact that nearby lenses have 
smaller values of v and r E owing to the smaller projection factor 
(see Section B3l . 

It has been suggested (e.g., Smith et al. 2002b) that parallax 
events may be preferentially caused by so-called disk-disk lens- 
ing, which refers to events in which both the source and lens lie in 
the disk. For such events the lens, source and observer are all co- 
rotating with the Galactic disk, producing small relative transverse 
velocities and hence enhancing the probability of detecting parallax 
signatures. This phenomenon can be seen in the right-hand panel 
of Fig. which exhibits a peak in the fraction of parallax events 
for nearby sources. However, the ratio of parallax events with disk 
sources is reduced due to the fact that the overall density of sources 
in the near side of the disk is only small. For our Galactic model 
we predict that ~5 per cent of parallax events have lensed sources 
lying within 5 kpc, which is a significant fraction when one consid- 
ers that only 0.5 per cent of all microlensing events have sources 
located in this region. 

In Table |4| we compare the relative abundance of parallax 



events within the different lens-source configurations, such as for 
disk-bulge events (i.e. events in which the lens is in the disk and 
the source is in the bulge), etc. Clearly, one can see that the major- 
ity of parallax events are due to the disk-disk and disk-bulge events, 
namely ~ 33 and ~ 38 per cent respectively. This can be contrasted 
with the most common configuration for all microlensing events, 
which is bulge-bulge. Note that our parallax events have ~ 5 per 
cent of sources within 5 kpc, which is slightly larger than the esti- 
mate of less than 3 per cent that was predicted by the simulations 
of Bennett et al. (2002a). 

4.3 Projected Einstein Radius and Projected Lens Velocity 

Fig.|S|plots the distribution of the Einstein radius projected into the 
observer plane r E and the projected velocity v for all events and 
for parallax events. This shows that parallax signatures are more 
readily detectable for events with small values of r E and v, as pre- 
dicted by Smith et al. (2002a). This can be understood, on consider- 
ing that the projected velocity of the lens should be comparable to 
(and preferably less than) the orbital velocity of the Earth, which is 
v e ~ 30 km s . The reason why smaller values of r E are favoured 
is because this parameter determines the length scale on which the 
magnification is calculated. This means that the magnitude of the 
deviations is determined by the magnitude of the Earth's motion 
relative to this projected Einstein radius, i.e. in general, the larger 
the projected Einstein radius, the smaller the deviations. In Fig.[9] 
we show the joint distribution of ? E and v, which - as expected - 
shows that parallax events preferentially have large / E and small v. 
From this figure, it is clear that the ability to detect parallax sig- 
natures is controlled by the quantity v/£ E = r E //|. For a given ac- 
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Lens 


Source 


Percentage of events with 


given lens-source configuration 






All events 


Weak parallax 


Moderate parallax 


Strong parallax 


near disk 


near disk 


3.1 (0.5) 


19.4 (3.1) 


33.2 (5.2) 


36.3 (5.7) 


near disk 


far disk 


12.0 (0.9) 


9.4 (3.4) 


7.5 (4.7) 


8.9(5.1) 


far disk 


far disk 


1.5 (0.1) 


0.9 (0.0) 


0.5 (0.0) 


0.8 (0.0) 


disk 


bulge 


29.6(17.1) 


39.6 (49.7) 


38.4 (62.8) 


36.9 (65.9) 


bulge 


bulge 


34.0 (71.6) 


14.3 (31.3) 


8.9 (17.4) 


7.3 (13.3) 


bulge 


disk 


19.9 (9.8) 


16.3 (12.4) 


11.6(9.9) 


9.8 (10.0) 



Table 4. The percentage of events as a function of lens-source configuration for all events and weak-, moderate-, and strong-parallax events. For the values 
without parentheses, a disk star is defined as one drawn from the population described by equation J4j, while a bulge star is drawn from equation J5}. Bennett 
et al. (2002a) used a different nomenclature, defining all stars within 3.5 kpc of the Galactic centre as "bulge" and all stars outside 3.5 kpc of the Galactic 
centre as "disk". To enable comparison with Bennett et al.'s results, we give in parentheses the values corresponding to these definitions. Note that the disk-disk 
events are subdivided according to whether the disk stars reside on the near- and far-side of the Galactic centre. 
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Figure 9. Contour plot showing the joint distribution of event time-scale 
and projected lens velocity on the observer plane for all events (dotted) and 
parallax events (solid). The four diagonal lines indicate constant v/t E = 
'W'e = '0.0, 1,0.1,0.01 km s day -1 . If we approximate the Earth's pro- 
jected acceleration to be o k 0.3 km s -1 day -1 (i.e. averaged over the whole 
year), then these contours correspond to J?l = a/(v/fE) ~ 0.03,0.3,3,30. 
This dimensionless acceleration parameter J?I effectively controls the sig- 
nificance of the parallax deviations (see Section l4~3l . 



celeration a, the dimensionless parameter = a/(rE/t E ) describes 
the deviation of the trajectory from a straight line (Smith, Mao & 
Paczyhski 2003b). Parallax deviations are caused by the Earth's or- 
bital acceleration. Although the magnitude of this acceleration is 
constant in the ecliptic plane, when projected into the lens plane 
it is no longer constant and so a varies throughout the year. How- 
ever, most parallax events last for a sufficiently long duration for 
the average value of a to be used. This means that the quantity 
r E /r| effectively determines the value of and hence controls the 
significance of the deviation due to the Earth's motion, i.e. a larger 
re/?! indicates that the parallax effect will be more difficult to de- 
tect and vice versa. We return to the issue of the f E distribution in 
the discussion, where we compare our predictions to the observed 
distribution. 



For a parallax event, it is possible to determine the orienta- 
tion of v, from which we can gain additional information regarding 
the event. In Fie. 1101 we show the distribution of this angle trans- 
formed into Galactic coordinates (9 a ), where 6 a is measured from 
the Galactic plane (in the direction of rotation) towards the North 
Galactic pole. From this figure, it can be seen that parallax events 
preferentially have v orientated in a direction parallel to the ecliptic 
plane. This can be understood when one considers that trajectories 
perpendicular to the ecliptic plane will have less time to be affected 
by the Earth's motion, compared to trajectories that pass parallel 
to the ecliptic plane, i.e. for perpendicular trajectories the Earth's 
acceleration will only affect the peak of the light curve (Bennett 
et al. 2002a). The orientation of v is also affected by the fact that 
parallax events are commonly due to bulge sources being lensed by 
foreground disk stars (see Section H^l . For the disk-bulge config- 
uration v is, on average, orientated along the direction of rotation 
of the Galactic plane owing to the fact that the Sun and the lenses 
share a common motion in the plane of the disk. This enhances the 
fraction of parallax events with -90° < 9 a < 90°, as can be seen 
fromFig.fTfjl 

We also investigated the radial velocities of the source stars for 
our mock catalogues. No significant difference was found between 
the parallax events and all lensed sources. 

4.4 Mass Estimators 

Bennett et al. (2002a) and Agol et al. (2002) proposed very similar 
techniques for estimating the mass of the lens for parallax events. 
Using our mock catalogues, we can assess the reliability of their 
estimator. 

Agol et al. (2002) assume that the source and lens popula- 
tions are both characterised by Gaussian velocity distributions with 
means (vl) and (vs) and dispersions 07, and as, diagonalised in 
Galactic longitude and latitude (£, b). Specialising to the case of a 
power-law mass function n(m) oc m~P, Agol et al.'s (2002) maxi- 
mum likelihood estimator becomes: 



L(x\v, t E ) cc /-»(1 - x) 5 -"^ exp 

0"fO"A 



20-2 . 

L,b / 



(14) 



where 

(y { , v b ) = (? L ) - *(vs) - (1 - x)(v Q + v), (15) 
and 

^ = oi/ + ^/. a i = < T l,b + xl < T l,b- (16) 
L(x\v, ?e) gives the likelihood of a lens lying at a fractional distance 
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Figure 10. The lower panel shows the distribution of the projected velocity angle, %, for all events (dotted line) and moderate parallax events (solid line) 
in our mock catalogue. 8a is measured in Galactic coordinates from the Galactic plane (in the direction of rotation) towards the North Galactic pole. The 
inset shows the projected velocity distribution in Galactic coordinates, where the horizontal axis corresponds to the direction of rotation and the vertical axis 
corresponds to the North Galactic Pole. The contours correspond to the 68.3% and 95.4% confidence intervals. The upper panels show the fraction of events 
that are classified as moderate parallax events as a function of the projected velocity angle. The dashed lines denote the orientation of the ecliptic plane in this 
coordinate system. 



x = D^/Ds, given the observables v and t^. This is obtained using 
Bayes' theorem, assuming uniform priors in lens distance and no 
errors on the velocity and timescale measurements. It can be con- 
verted to a likelihood in mass m using 



m(x) 



v 2 ?|c 2 (l - x) 



(17) 



As the mass functions in our simulations are more complicated than 
simple power-laws, we explore two different choices, namely /? = 
and /i = 1.5. In the latter case, equation 1 141 reduces to a formula 
proposed by Bennett et al. (2002a). 

The estimator also assumes that the distance of the source D s 
is known and that the lensing population is identifiable, so that an 
informed choice for the density p L can be made. We extracted 728 
events from our mock catalogue that pass Wozniak et al.'s (2001) 
criteria. These are moderate parallax events with A^ 2 > 50 and 
S IN > 5. We use the fitted values of v, /e, together with the known 
source distance D s and the known type of deflector population to 
compute the estimated mass. Of course, the accuracy in the estimate 
can be calculated, as the true mass is known. 

Fig. II II shows histograms of the logarithm of the ratio of the 
recovered mass to the true mass for the 89 events caused by bulge 
deflectors (left panels) and 639 events caused by disk deflectors 
(right panels) using the likelihood estimator with yS = 1.5 and ft = 0. 
For the bulge lenses, both estimators give a similar performance 
with mean percentage errors of -22 per cent (ft = 0) and -45 per 
cent (p = 1.5). A negative value of the mean percentage error im- 
plies that the estimators typically underestimate the mass. For the 



disk lenses the estimators tend to overestimate the mass, with a 
mean percentage error of 85 per cent for f} = 0, compared to 20 per 
cent for fj = 1.5. The simulations use more complicated mass func- 
tions than power-laws. It is noteworthy that the effect of using a 
simple power-law mass function in the likelihood estimator (which 
is needed to perform the integrals) is to cause systematic offsets in 
the recovered masses as compared to the true masses. 

The interpretation of the effectiveness of these mass estimators 
suffers from an additional complication. As was shown in Section 
13.11 the accurate recovery of the parallax parameters such as v is 
not guaranteed, even for our convincing parallax events (see Fig. 
[4}- Clearly these mass estimators cannot be expected to perform 
reliably if the parallax parameters have not been determined accu- 
rately from the observed light curves. 



4.5 Long Duration Events and Black Hole Candidate Lenses 

Our simulations can also be used to investigate the nature of ob- 
served long duration microlensing events, a number of which have 
been speculated to be caused by black hole lenses. 

In our model, we do not include any stellar remnant pop- 
ulations. However, white dwarfs (M ~ 0.6 Mq), neutron stars 
(M ~ 1.4 Mq) and black holes (M ~ several Mq) must exist in the 
Galactic disk, and so their contribution to parallax events is not ac- 
counted for in our model. For long durations (t£ <; 100 d), the time- 
scale distribution asymptotically approaches a power-law (Mao & 
Paczyhski 1996). As a result, the fraction of events contributed by 
stellar remnants reaches an asymptotic value, determined entirely 
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Figure 11. The error distributions of the estimators of the mass of the lens of a parallax event using /? = 1 .5 (above) and /J = (below) in equation 1141 . The 
results for bulge deflectors are shown on the left, those for disk deflectors on the right. The abscissa is the logarithm of the ratio of the recovered to the true 
mass. 



by the mass function. The relative fraction is weighted according 
to ~ M 2 n(M)dM, and hence favours massive lenses (Agol et al. 
2002). If we adopt a mass function for stellar remnants (e.g., Gould 
2000; Han & Gould 2003), we find that the fraction of long du- 
ration events contributed by stellar remnants is about 56 per cent. 
Therefore our predicted parallax fraction may be too low by a fac- 
tor of ~2. However, in reality, the discrepancy may be less than this 
because on average massive stellar remnants will have larger f E and 
correspondingly a reduced probability of exhibiting parallax signa- 
tures (see Fig. [8] and Section |4~3*1 . In addition, since the baseline 
of our mock light curves is only 3 yr, events with such long time- 
scales may not have been considered in this analysis. For a mock 
event to be included in our parallax catalogue, we require a constant 
baseline, which excludes exceptionally long duration events. 

To date, the most promising candidate black hole event is 
OGLE-1999-BUL-32/MACHO-99-BLG-22 (Mao et al. 2002; Ben- 
nett et al. 2002b). This event, along with an additional two 
events MACHO-96-BLG-5 and MACHO-98-BLG-6 (Bennett et al. 
2002a), have been analysed by Agol et al. (2002). In this paper, they 
used their mass estimator (see equation^} to conclude that event 
OGLE-1999-BUL-32 has a probability of 76 per cent of being 
caused by a black hole lens (with the remaining two events having 
probabilities less than 20 per cent). The event parameters t E and f E 
for these three events can be found in TableQ In addition to t E and 
r E , Agol et al. (2002) also incorporated the angle of the projected 
velocity 6 a and the lower limit on the lens magnitude into their es- 
timator. For events OGLE-1999-BUL-32, MACHO-96-BLG-5 and 
MACHO-98-BLG-6, the respective angles and /-band lens magni- 
tude constraints are: 142°, I L > 18.6; -20°, I L > 18.6; -68°. The 
lower limit for the /-band lens magnitude of event MACHO-98- 



BLG-6 is unknown since Bennett et al. (2002a) only constrain the 
V-band magnitude. 

Using our simulations we now attempt to determine whether 
the observed characteristics of these events are consistent with our 
stellar population (that contains no remnants). In Fie. 1121 we cal- 
culate the joint probability distribution of the parameters f E and t E 
for the three black hole candidate events; for each event, we plot 
the probability distributions for the mock events that match the ob- 
served do and have lens magnitudes fainter than the required limit. 
Note that here we do not incorporate the OGLE-II event detection 
efficiency, since many of these long duration events will be omit- 
ted due to the 3-yr baseline of the experiment. From this figure, it 
can be concluded that although the properties of event MACHO- 
1998-BLG-6 are consistent with what is predicted from our sim- 
ulations, the remaining two events, especially OGLE-1999-BUL- 
32, are clearly inconsistent. However, such large values of f E can 
be obtained by having the lens and source close together, such as 
from bulge self-lensing. Although bulge self-lensing events typi- 
cally have a shorter time-scale (with a mean t E of 26.3 d versus 32. 1 
d for all events), if the Galactic model incorporated the streaming 
motion of the bar then it may be feasible to obtain both large f E and 
t E simultaneously. 



5 DISCUSSION AND CONCLUSIONS 

This paper has investigated the properties of parallax microlens- 
ing events towards the Galactic bulge. The expected fraction of 
parallax events based on standard models of the Galaxy is of the 
order of a few per cent, which is compatible with most observed 
findings. However, it should be noted that observed catalogues fea- 
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Figure 12. Contour plot showing the joint distribution of event time-scale and projected Einstein radius on the observer plane for all events regardless of 
whether or not they pass the OGLE-II detection criteria of Wozniak et al. (2001). The contours denote the 68.3, 95.4 and 99.7 per cent confidence regions. 
Each panel corresponds to the mock events that have projected velocity angle 0q and lens magnitude consistent with the observed black hole candidate named. 
Note that no /-band magnitude limit is available for event MACHO-1998-BLG-6 and hence we provide two sets of contours: the solid line shows the effect 
of including no restriction on the lens magnitude and the dotted line corresponds to the limit /l > 21.2, which is obtained from the V-band limit VL > 22.2 
(Bennett et al. 2002a) and a fiducial estimate of the colour of the lens (V - /)l = 1.0. Note that since one might naively suspect that (V - /)l <: 1.0, the true 
contours should lie somewhere between the dotted and solid contours presented here. 



ture other types of exotic microlensing events, such as binary-lens 
events (see Jaroszyhski 2002) and binary-source (xallarap) events, 
many of which can mimic parallax signatures and hence display an 
improvement in x 2 when fit with a parallax model (e.g. Smith et 
al. 2002b). Although our predicted parallax fractions seem consis- 
tent with observed catalogues, there is one exception, namely the 
results from the EROS collaboration (Afonso et al. 2003). In this 
paper, they report a fraction approximately 12.5 per cent, which 
is even more surprising when one considers that this catalogue is 
from clump-giant sources, which supposedly reside in the bulge; 
our simulations predict that the fraction of parallax events for bulge 
sources is only ~ 1 per cent. It could be argued that the EROS 
source stars are bright red clump giant stars and as such may have 
higher S/N than our simulations. However, even if we restrict our 
simulations to bright bulge source stars with high S/N, the percent- 
age of events with convincing parallax signatures is only 2 per cent; 
given this predicted fraction, the probability that such a sample of 
16 events will yield two parallax events is only 0.04. Therefore, 
although no strong conclusions can be made, it seems that this ob- 
served fraction of parallax events from the EROS catalogue may be 
inconsistent with our simulations. 



Our models have one major shortcoming, namely that we only 
included blending from the lenses themselves, i.e. we did not in- 
clude blending from nearby unrelated stars. We attempted to limit 
this problem by only considering sources brighter than / = 19, 
since it is commonly assumed that such bright sources are less af- 



fected by blending. 4 To test what effect blending from unrelated 
stars would have on our simulations, we perform the following 
simple investigation; we generate an additional catalogue of mock 
events with blending ratio (i.e. ratio of source flux to total baseline 
flux) distributed uniformly between zero and one. This means that 
the source stars for the blended catalogue are fainter than the cor- 
responding source star in the unblended catalogue. Although the 
mean time-scale of observed events (t E ) is practically unchanged in 
this new catalogue, the overall fraction of parallax events is slightly 
reduced due to the greater number of lower signal-to-noise ratio 
events, which can be expected since the source stars are fainter for 
the blended catalogue. Currently, it is unclear how blending varies 
as a function of source magnitude and event time-scales, so we 
cannot address this question quantitatively. However, the qualita- 
tive differences between all microlensing and parallax events (i.e. 
the distributions shown in Section^) appear unchanged from this 
analysis of blending. 

To obtain a greater understanding of the underlying nature of 
parallax events, we examined the distributions of the event param- 
eters r E and v (or, equivalently, r E ) for parallax events, comparing 
them to the distributions for all events. One important question is 
how well do our simulations match the observed characteristics of 
parallax events. Although the overall fractions of parallax events 
seem consistent with current catalogues, we can also tentatively 
compare the parallax parameters such as the projected Einstein ra- 
dius, f E . In Fig.^| we show a plot comparing the predictions from 

4 However, recent work by Sumi et al. (2005) has suggested that even this 
assumption may not be wholly reliable. 
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Figure 13. A comparison between the distribution of the projected Einstein 
radius (f E ) for observed parallax events and our mock catalogues. The dis- 
tributions from our mock catalogue are given by the solid lines and the ob- 
served distributions are given by the data points (where Poisson errors have 
been assumed). The lower panel shows the observed distribution, while the 
middle panel shows the underlying distribution, after correcting for the ef- 
ficiency (top panel). 



our simulations with the distribution of f E for the observed paral- 
lax events from Table Q It should be noted that since the events 
from TableQcome from a wide range of microlensing experiments 
with differing durations, sampling and photometric properties, any 
comparison can only be very speculative. However, as can be seen 
from the bottom panel of this figure, the observed distribution of 
f E (as given by the data points) appears to be in reasonable agree- 
ment with the distribution from parallax events in our mock cat- 
alogues. Furthermore, from our simulations, we can calculate the 
efficiency of recovering parallax signatures and attempt to deter- 
mine the form of the underlying r E distribution for these observed 
events. This method is analogous to the one used to convert an ob- 
served time-scale distribution into a true underlying time-scale dis- 
tribution using the 'detection efficiency'. The upper panel of Fig. 
IQI shows the efficiency of detecting parallax signatures, and when 
applied to the observed f E distribution we obtain the underlying 
distribution shown in the middle panel. It appears that the observed 
distribution is consistent with the underlying f E distribution from 
our simulations. 

Although there are currently too few observed parallax events 
for us to make any firm statements about the distribution of the 
parallax parameters, it may become possible in the near future for 
projects such as OGLE-III (Udalski 2003), which is currently de- 
tecting ~500 microlensing events each year. Assuming the rate of 
parallax detection for this project is similar to that of OGLE-II, the 
three-fold increase in sky coverage means that a few good qual- 
ity parallax events will be detected each year. It is expected that 
the OGLE-III project (or an upgraded experiment, OGLE-IV) will 
continue for many years, which should be sufficient to enable the 
construction of distributions of r E and v that are reasonably well- 
constrained. A similar study to the present one would allow a de- 



tailed comparison between theoretical predictions and observations 
for parallax events. 

Unique lens mass determinations for parallax events may be- 
come routine in the near future by combining a reliable measure- 
ment of r E with additional constraints, such as has been shown al- 
ready for a number of events (An et al. 2002; Gould et al. 2004; 
Kubas et al. 2005). One particular approach that has been pro- 
posed is the measurement of the separation between the two mi- 
crolensed images using powerful interferometers such as VLTI 
(Delplancke, Gorski & Richichi 2001). We can see from our sim- 
ulations that such image separations at maximum magnification 

(#se P = S E yjug + 4, where u Q is the impact parameter corresponding 
to the maximum magnification) are larger for our parallax sample, 
with mean separations of 1.8 mas compared to 0.8 mas for the full 
sample. Therefore it will be easier to resolve the two microlensed 
images for parallax events and, in addition, their longer duration 
will make it more feasible to plan high signal-to-noise ratio obser- 
vations while the events are still undergoing high magnification. 
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